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基于 GF-1/WFV 时 间 序 列 的 葡萄 识别 模型 
-以 宁夏 红寺堡 区 为 例 " 
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(1. 宁夏 大 学 西北 土地 退化 与 生态 系统 恢复 省 部 
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基建 国家 重点 实验 室 培 育 基地 ,宁夏 银川 750021; 


2. 宁夏 大 学 西北 退化 生态 系统 恢复 与 重建 教育 部 重点 实验 室 ,宁夏 银川 750021; 


3. 南京 大 学 国际 地 球 系统 科学 研究 所 ,江苏 ”南京 ”210093) 


摘 X: 以 宁夏 红寺堡 区 为 研究 区 ,基于 高 分 一 号 (GF - 1XWFV ) 卫星 构建 葡萄 生长 季 时 间 序 列 光 谱 数 据 , 运 用 


(Jeffreys-Matusita)(J-M) 距 离 分 析 和 葡萄 地 块 归 一 化 植被 指数 (NDVI) 时 序曲 线 特 征 确定 了 最 佳 识别 时 相 , 将 最 佳 时 


TH 


相 的 NDVI、 相 邻 时 相差 值 速率 和 曲线 积分 训练 样本 集 导 入 Clementine 数据 挖掘 软件 中 ,利用 C5. 0 决策 树 分 类 算 


法 ,并 结合 专家 经 验 法 构建 葡萄 林 决 策 树 提取 模型 。 


结果 表明 :构建 的 识别 模型 能 够 满足 葡萄 的 识别 需求 ,但 在 不 


同 履 盖 度 的 葡萄 地 块 上 精度 有 所 差异 ;基于 决策 树 分 类 的 总 体 精 度 为 93.7196 ,Kappa 系数 为 0.91 。 其 中 ,中 低 履 
盖 度 葡萄 林 生 产 精度 为 90. 82% ,用 户 精度 为 88.56% ;高 覆盖 度 葡 萄 林 生 产 精度 为 92. 44% ,用 户 精度 为 91. 18% 。 
关键 词 : 葡萄 林 ; 遥感 提取 ; GF - L/WEV 时 序数 据 ; 识别 模型 ;曲线 积分 ; 决策 树 ; 宁夏 


简 萄 是 分 布 最 为 广泛 的 果树 树种 之 一 ,如 何 快 
速 获取 葡萄 的 空间 分 布 及 其 生长 状态 ,为 葡萄 种 植 
的 精细 化 管理 .产业 的 可 持续 发 展 及 其 生态 效应 评 
价 提供 技术 支撑 ,是 当前 急需 解决 的 问题 。 农 业 通 
感 技术 的 发 展 使 得 快速 ,准确 监测 农作物 空间 分 布 
得 以 实现 ,对 于 农业 的 生产 与 布局 调控 具有 重大 意 
义 " ,可 作为 开展 葡 菊 种植 时 空 分 布 调查 的 重要 手 
段 。 宁 夏 资 兰 山 东 荔 和 葡萄酒 产 区 是 我 国 八 大 葡萄酒 
产 区 之 一 ,其 中 ,红寺堡 区 是 种 植 优质 酿酒 葡萄 的 最 
佳 产 区 之 一 ”"。 红 寺 堡 区 作为 新 绿洲 开发 区 , 摸 清 
葡萄 种 植 区 的 具体 分 布 情况 对 于 当地 葡萄 的 合理 种 
植 及 管理 极为 重要 。 

内 外 学 者 利用 遥感 技术 在 农作物 信息 的 识别 
和 提取 方面 已 进行 了 大 量 研究 工作 。 目 前 ,用 于 作 
物 面积 监测 的 遥感 数据 源 以 多 光谱 影像 为 主 , 按 空 
间 分 辨 率 的 不 同 分 为 高 .中 高 和 中 低 分 辩 率 。 高 空 
间 分 辨 率 影像 虽然 精度 高 ,但 时 间 分 辨 率 过 低 ,不 能 
准确 获取 满足 应 用 需求 的 作物 关键 生育 期 影像 ; 低 
空间 分 辨 率 的 遥感 影像 虽 满 足 时 间 序 列 的 实践 需 
求 , 但 空间 分 辩 率 非常 差 , 只 能 在 一 定 精度 范围 内 进 
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行 耕 地 提取 ”) ;而 中 分 辩 率 可 以 弥补 高 空间 分 辩 率 
和 低空 间 分 辩 率 的 缺陷 ,并 广泛 应 用 于 作物 面积 提 
取 以 及 耕地 提取 等 "*“。 虽 然 中 分 辨 率 影 像 中 混合 
像 元 现象 较 多 ,地 物 识别 性 较 弱 ,但 充分 利用 计算 机 
分 类 方法 ,辅助 于 地 物 在 光谱 中 呈现 的 规律 可 提高 
提取 精度 。 其 中 ,基于 数据 挖掘 技术 的 决策 树 分 
类 法 可 以 处 理 多 维 属性 时 确定 各 个 属性 的 重要 程 
BE ,提取 重要 性 高 的 属性 作为 分 类 判 据 并 生成 分 类 
规则 ,在 遥感 影像 的 分 类 和 目标 识别 当中 应 用 较为 
pue 

为 了 更 准确 提取 作物 种 植 信息 , 现 已 逐渐 由 利 
用 单一 时 相 影 像 向 利用 多 时 相 影 像 的 模式 转变 ,多 
时 相 影 像 能 够 有 效 监 测 不 同 作物 生长 情况 ,更 利于 
作物 的 区 分 。 因 此 ,时 间 序 列 的 构建 逐渐 成 为 热 
m, HÆK EK 小麦 .棉花 等 农作物 提取 方面 研 
究 广泛 ”” ,但 在 果林 提取 方面 相对 较 少 ,在 葡萄 
林 提 取 方 面 的 研究 更 鲜 见 报道 。 我 国正 在 实施 的 
“高 分 辩 率 对 地 观测 ”系统 重大 专项 中 ,高 分 辩 率 系 
列 卫 星 中 的 首 颗 光 学 卫星 高 分 一 号 (CF -1) 已 在 轨 
运行 近 2 a, 获 取 了 大 量 高 质量 数据 , 正 处 于 应 用 人 研 
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究 阶 段 ,其 能 够 在 自然 资源 调查 、 植 被 变化 监测 、 环 
境 质量 评估 及 生态 环境 变化 等 研究 中 发 挥 重要 作 
用 -中 。 此 外 ,GF -1 卫星 具有 高 .中 空间 分 辩 率 
对 地 观测 和 大 幅 宽 成 像 结合 的 特点 ,同时 搭载 了 幅 
宽 优 于 60 km 的 2 m 分 辩 率 全 色 和 8 m 分 辨 率 多 光 
谱 传 感 器 , 幅 宽 优 于 800 km 的 4 & 16 km 分 辨 率 多 
光谱 传感器 ,理论 上 可 以 实现 每 4 d 提供 1 次 数据 ， 
构成 时 间 序 列 数据 ,从 而 实现 高 空间 分 辨 率 和 高 时 
间 分 辩 率 相 结合 的 独特 优势 ,这 也 是 当前 遥感 对 地 
观测 技术 发 展 的 重要 方向 "9 。 因 此 ,开展 基于 
GF -1 时 序数 据 的 宁夏 葡萄 主 产 区 种 植 空 间 分 布 ， 
对 于 遥感 典型 植被 类 型 识别 与 监测 研究 领域 新 技术 
的 研发 具有 重要 意义 。 

本 文选 择 兼 具 较 好 时 间 和 空间 分 辨 率 的 中 高 分 
辨 率 卫星 高 分 一 号 ,作为 葡萄 林 基 于 时 间 序 列 提取 
提供 较 好 的 数据 源 ,结合 遥感 影像 的 光谱 特性 和 提 
取 的 植被 指数 运用 决策 树 分 类 算法 ,提取 和 识别 葡 
萄 林 空 间 分 布 信息 。 以 宁夏 贺兰山 东 苑 葡萄 酒 产 区 
的 红寺堡 区 为 研究 区 ,基于 GF -1AWFV 时 序数 据 ， 
采用 决策 树 分 类 方法 ,分 为 中 低 履 盖 度 葡萄 林 和 高 
覆盖 度 葡 萄 林 , 分 别提 取 宁 夏 红 寺 堡 区 葡萄 林地 , 探 
索 并 建立 普 适 性 与 操作 性 较 强 的 葡萄 林地 决策 树 分 
类 模型 ,从 而 为 葡萄 林地 提取 提供 参考 ,也 为 未 来 葡 
萄 种 植 合理 布局 及 其 生态 环境 保护 提供 借鉴 。 


1 数据 与 方法 


1.1 研究 区 概况 

红寺堡 区 位 于 贺兰山 东 荔 葡萄酒 产 区 的 最 南 
端 ,地 理 坐 标 为 105945' ~ 106° 30' E, 37? 10' ~ 
37?29'N ,地 势 由 东南 向 西北 倾斜 。 该 区 域 属 典 型 温 
带 大 陆 性 气候 ,干旱 少雨 ,多 年 平均 降水 量 为 251 
mm, 由 东南 向 西北 递减 ,降雨 多 集中 在 6 一 9 月 , 占 
全 年 降水 量 的 72. 4% 。 多 年 平均 蒸发 量 为 2 387 


mm, im H £23 13.7 *C ,全 年 日 照 时 数 为 2 900 ~ 
3 550 h。 该 区 土地 总 面积 约 为 880.73 km ,良好 的 
气候 条 件 极 有 利于 葡萄 种 植 ,近年 来 ,葡萄 产业 已 成 
为 红寺堡 的 第 一 支柱 产业 沁 - 。 因 此 ,以 红寺堡 区 
的 葡萄 林 为 研究 对 象 进行 葡萄 林 识 别 研究 极 具 代表 
性 和 典型 性 。 
1.2 数据 源 与 预 处 理 

高 分 一 号 (GF - 1) 卫星 搭载 的 4 台 多 光谱 相机 
(WFV) 幅 宽 可 达 800 km, 包 括 蓝 (0.45 ~ 0. 52 
um) 绿 (0.52 ~0.59 hm)、 红 (0.63 ~0.69 um), 
近 红 外 (0.77 ~0. 89 pm)4 个 波段 ,空间 分 辨 率 达 
到 16 m , 重 访 周期 为 2 d, 同 时 具有 较 高 空间 和 时 间 
分 辨 率 的 特点 ,能 获取 无 云 的 植被 关键 生育 期 的 遥 
感 影像 。 考 虑 到 所 选 影像 应 该 在 同一 天 内 履 盖 红 寺 
EX, 且 影 像 应 该 无 云 或 少 云 , 共 选 取 了 2017 年 
3—12 月 共 12 个 时 相 的 遥感 数据 ( 表 1) 。 

影像 的 预 处 理 包 括 辐 射 定 标 ` 大气 校正 和 几何 
校正 ,其 中 大 气 校正 采用 Flaash 模型 。 最 后 裁剪 拼 
接 最 终 获 得 GF -1/WFV 的 时 间 序 列 数据 集 。 利 用 
多 时 相 NDVI 数据 能 够 充分 刻画 出 同一 作物 在 不 同 
生育 阶段 的 光谱 差异 ,可 为 准确 提取 作物 种 植 面积 
信息 提供 技术 支持 ""。 研 究 区 时 序 归 一 化 植被 指 
Ži ( normalized difference vegetation index, NDVI) 计 
算 公 式 如 下 : 


NDVI = ( Nir - Red )/ ( Nir + Red) (1) 
式 中 : Nir 为 近 红 外 波段 反射 率 ; Red 为 红 波 反射 
率 。 由 于 本 研究 采用 的 GF -1 数据 为 还 感 影像 , 故 


近 红 外 波段 为 Band 4, 红 光波 段 为 Band 3。 

为 获取 宁夏 红寺堡 灌区 主要 地 物 类 型 的 分 布 ， 
2017 年 8 月 对 红寺堡 区 进行 地 面 实 地 调查 (图 1)。 
野外 调查 时 采用 手持 GPS ( global positioning system) 
测量 特征 地 块 的 经 纬度 坐标 并 记录 植被 类 型 。 此 次 
调查 共 获 取样 本 点 1 901 个 ( 表 2) 。 


表 1 影像 时 间 序 列 信息 
Tab.1 Time series information of images 
序号 传感器 类 型 和 产品 序号 影像 获取 时 间 | 序号 传感器 类 型 和 产品 序号 影像 获取 时 间 
1 GF1_WFV4_E106.7_N38.5_L1 A0002263952 2017 -03 -26 7 GFI WFV3 E106.1 N38.9 L1A0002545275 2017 -08 -16 
2 GFI1 WFV4 E106.7 N36.7 L1A0002331149 A 2017 -04 -28 8 GFI WFV4 E106.6 N38.5 L1A0002580944 2017 -09 -02 
3 GF1 WFV3 E106.2 N38.9 L1A0002379084 . 2017 -05 -26 9 GFI WFV2 E107.5. N39.3 L1A0002671474 2017 -10 -12 
4 GFI WFV4 E106.7 N38.5 L1A0002424156 2017 -06 -16 10 GFI WFV2 E107.1. N39.3 L1A0002755542 2017-11-10 
5 GF1 WFV3 E106.1 N38.9 L1A0002464407 . 2017 -07 -06 11 GF1_WFV3_E107.2_N38.9_L1A0002767349 2017-11-14 
6 GFI WFV1 E105.5 N38.0 L1A0002515969 . 2017 -07 -30 12 GFI WFV3 E106.4 N38.9 L1A0002824873 2017 -12 -05 
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图 1 宁夏 红寺堡 位 置 及 样 点 分 布 示 意图 


Fig.1 Schematic diagram of the location of Hongsibao and the distribution of sample sites 


表 2 地 类 取样 点 个 数 


Tab.2 Number of sampling sites of land use types 


地 类 训练 样本 点 数 ”验证 样本 点 数 ”样本 总 数 
中 低 覆 盖 度 葡萄 林 215 155 370 
高 覆盖 度 葡萄 林 135 92 227 
玉米 181 100 281 
小 麦 117 95 212 
牧 革 和 蔬菜 96 70 166 
防护 林 106 65 171 
fii C Hl 137 84 221 
裸 地 137 120 257 


2 研究 方法 与 技术 路 线 


CL) 根据 研究 区 葡萄 生长 物候 历 , 确 定 葡萄 生 
长 的 关键 时 期 ,选取 相应 时 间 影 像 ;(2) 影像 在 图 像 
预 处 理 后 ,利用 ENVIS. 3 软件 计算 所 选 11 个 时 相 
的 NDVI 值 ,构建 完整 时 间 序 列 的 GE/WFV - NDVI 
数据 ,通过 分 析 得 到 葡萄 提取 的 最 佳 时 序 组 合 ; (3 ) 
将 NDVI 数据 进行 波段 运算 ,分 别 计算 其 相 邻 时 相 
差 值 速率 (K) 和 曲线 积分 (S) ;(4) 将 NDVI、 相 邻 时 
相差 值 速率 (K) 和 曲线 积分 (5) 训练 样本 集 导 入 
Clementine 数据 挖掘 软件 中 ,利用 C5.0 决 策 树 分 类 
算法 ,并 结合 专家 经 验 法 构建 葡萄 林 决 策 树 提 取 模 
型 ,进行 研究 区 葡萄 种 植 区 域 识别 ,并 进行 精度 验 
证 。 本 文 的 研究 技术 路 线 如 图 2 所 示 。 
2.1 葡萄 生长 周期 

研究 区 冬季 寒冷 ,为 预防 葡萄 茧 风干 和 受 冻 , 需 
对 其 进行 埋 土 防寒 "”] 。 翌 年 春季 4 月 中 旬 将 其 挖 
出 ,5 月 初 开始 萌芽 抽 枝 ,5 月 20 日 左右 进入 花期 ;6 
月 中 旬 ( 打 春 梢 ) 左 右 进入 坐果 期 ;7 一 8 月 葡萄 枝 蔓 
和 果实 进入 生长 旺 期 ,8 月 中 旬 在 枝 划 上 打 秋 梢 以 
减少 养分 流失 ; 9 一 10 月 葡萄 果实 成 熟 , 之 后 枝 草 开 


多 时 相 GF-1 WFV 数 据 
生成 数据 集 


NDVI 时 间 序 列 影像 集 训练 样本 集 
i] s 
典型 地 类 NDVI 
时 序曲 线 分 析 


到 2 技术 路 线 


Fig.2 Technology route 


RU XC ,采摘 完成 时 将 叶子 清除 ,只 留 主 葵 并 将 其 掩 
Hg ,其 后 葡萄 进入 休 眼 期 ( 表 3)。 
2.2 ”典型 地 物 NDVI 时 间 序 列 构建 

时 间 序 列 植被 指数 均值 变化 能 直观 反映 一 年 中 
同一 植被 在 不 同时 间 和 同一 时 间 不 同 植被 的 植被 指 
数 差异 。 归 一 化 植被 指数 与 植被 的 许多 参数 密 
切 相 关 , 如 叶 面 积 指 数 、 植 被 覆盖 度 和 植被 郁 闭 度 
等 ,利用 NDVI 能 够 很 好 地 提取 植被 并 且 具 有 较 强 
Fb THESE 77 7, 

由 于 葡萄 栽培 在 秋季 清理 枝 葛 和 冬季 掩埋 的 特 
殊 性 ,所 以 研究 中 选取 包括 秋冬 季节 的 全 年 时 间 序 
列 为 研究 周期 ,以 充分 展现 其 年 内 的 生长 变化 特征 ， 
探讨 其 提取 方法 。 人 研究 的 主要 工作 是 将 葡萄 林 与 耕 
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I 模型 一 一 以 宁夏 红寺堡 区 为 例 


RI 不 同时 相 遥 感 影像 与 葡萄 生长 周期 对 照 


Tab.3 Compared results of different-phase remote sensing images and grape growth cycle 


E ie 幼 果 期 à TE RAMO EU E EE 
生长 周期 04 -20— 05 -21— 06 - 16— 07 -25— 09 -01— 09 -26— 10 -21— 11 -25— 
H-H 05 -20 06 - 15 07 -20 08 -30 09 -25 10 -20 11 -20 04 -10 
GF-1 2017-04-28 2017 -05 -26 Ma ru pa 2017 -09 -02 2017 -10-12 2017-11-14 E 2d 
地 、 林 地、 裸 地 分 开 。 基 于 红寺堡 区 主要 的 地 物 类 NDVI 差 值 曲线 积分 . 
型 ,利用 公式 (1) 计 算 影 像 的 NDVI 值 , 求 取 每 种 
利用 公式 (1) 计算 影像 的 值 , 求 取 每 种 地 S, -0444 yx 2 T)2 (5) 


类 在 不 同时 间 的 NDVI 均值 ,建立 研究 区 主要 地 物 
NDVI 均值 时 间 序 列 曲线 。 
2.3 Jeffreys-Matusita 距离 计算 与 最 佳 时 相 

NDVI 时 间 序 列 共 有 12 个 时 相 , 包 括 葡 萄 掩埋 
期 和 生长 发 育 期 的 数据 ,其 数据 量 较 大 ,因此 ,可 进 
行 最 佳 时 序 组 合 分 析 。 研 究 的 最 终 目的 是 将 葡萄 林 
和 其 他 地 类 区 分 ,地 表 特 征 可 分 性 的 判定 方法 很 多 ， 
如 J-M 距离 (Jeffreys-Matusita distance ) , B. 距离 
( Bhatta-charyya distance ) 、 离 散 度 、 样 本 间 平 均 距 
离 ,类别 间 相对 距离 等 。 相 对 于 其 他 指标 ,J-M 距离 
被 认为 更 适合 表达 类 别 可 分 性 。 基 于 某 一 特征 两 类 
样本 的 J-M 距离 ,计算 公式 表示 为 : 

J (2) 

RPJ 表示 J-M 距离 ;。 表示 某 类 特征 的 均值 ;有 表 
示 在 某 一 特征 维 的 巴 氏 距离 。 在 样本 对 象 满足 正 态 
分 布 时 ,两 类 样本 对 象 的 巴 氏 距离 为 : 


5,45; 
26,6, 


1 2 1 
B-3g(6 Te) + + 2 n 


| (3) 


式 中 :6 表示 某 类 特征 的 方差 。 

J-M 距离 值 在 0 ~2 之 间 , 其 大 小 代表 样本 间 可 
分 离 程 度 的 高 度 。J-M 距离 大 于 1.8 为 合格 ,说 明 
样本 间 具 有 很 好 的 光谱 可 分 性 ,其 值 越 大 表明 地 物 
间 可 分 离 性 较 好 "1。 
2.4 NDVI 差 值 速率 与 相 邻 曲线 积分 计算 

季节 变化 作为 植被 覆盖 的 特征 之 一 ,也 是 不 同 
植被 类 型 随时 间 变 化 的 一 种 自然 现象 ,利用 不 同时 
相 条 件 下 的 NDVI 序列 ,可 以 比较 准确 地 反映 出 植 
物 生长 的 季 相 变化 。 在 NDVI BIER di AR PT, 
加 入 特殊 时 相 的 差 值 速率 及 曲线 积分 ,可 作为 进 一 
步 识别 相似 地 类 的 依据 。 

计算 公式 : 

A-A; 
NDVI 差 值 速率 :K =p y 


+ 


(4) 


式 中 :4 分 别 表示 各 时 相 的 NDVI; i 表示 NDVI 时 相 
序列 编号 ;7 为 相 邻 2 个 时 相 NDVI 的 天 数 差 ( 表 
4) ;K 表示 各 相 邻 时 相 的 NDVI 差 值 速率 ;S; 表示 各 
相 邻 时 相 的 NDVI 曲线 积分 值 。 


R4 NDVI 数据 及 编号 
Tab.4 NDVI data and serial numbers 


NDVI 编 号 (i) ”NDVI 时 相 序 列 (8) ”NDVI 对 应 天 数 (7)/d 
1 03 -26 85 
2 04 -28 118 
3 05 -26 146 
4 06 -16 167 
5 07 -06 187 
6 07 -30 211 
7 08 -16 228 
8 09 -02 245 
9 10 -12 285 
10 11 -14 318 
11 


12 -05 339 


2.5 基于 决策 树 的 葡萄 林 识 别 模型 

C5.0 作为 一 种 决策 树 算 法 ,能 根据 最 大 信息 增 
益 (information gain) ) 自动 计算 出 达到 最 佳 决策 树 
时 各 叶子 国 值 ,从 而 提高 决策 树 分 类 算法 的 客观 性 
与 精度 ”-” 。 因 此 ,本 文采 用 数据 挖掘 软件 Clem- 
entine 对 数据 集 深 入 挖 据 ,并 结合 专家 经 验 法 在 EN- 
VI5.3 中 生成 决策 树 模型 及 分 类 规则 。 

本 研究 主要 对 易 混淆 的 高 履 盖 植被 类 (防护 
PR HEHE V E83 6E i: SE 8] 28] ) 和 低 覆 盖 植 被 类 ( HP ICE 
盖 度 葡萄 林 jfi Ej). 及 裸 地 进行 遥感 识别 ,从 而 实 
现 葡 萄 林地 的 高 精度 提取 。 因 此 ,选取 人 研究 区 不 同 
植被 覆盖 地 类 ,分 别 分 析 NDVI 曲线 特征 及 其 差 值 
速率 特征 、 积 分 值 等 特征 变量 的 平均 值 , 找 出 葡萄 林 
与 其 他 不 同 植被 覆盖 地 类 所 存在 的 差异 性 ,从 而 实 
现 目标 地 类 提取 目的 。 


3 结果 与 分 析 


3.1 典型 地 物 GF - 1/WFV-NDVI 时 间 序 列 

从 图 3 可 以 看 出 ,中 低 覆 盖 度 葡萄 林 的 NDVI 
曲线 自 萌芽 期 至 花期 呈 陡 坡 形式 增长 ,并 于 果实 膨 
大 及 着 色 期 达到 最 大 峰值 , 受 8 月 中 旬 打 秋 梢 活动 
的 影响 , 随 生物 量 的 减少 ,NDVI 曲线 在 该 时 期 出 现 
低谷 ,至 果实 成 熟 期 后 葡萄 林 生 长 较为 旺盛 , NDVI 
曲线 达到 另 一 个 峰值 。 中 低 覆 盖 度 葡萄 林 较 明显 的 
曲线 特征 在 于 8 一 9 月 低 于 耕地 和 林地 , 出现 区 分 中 
低 覆 盖 度 葡萄 林 和 耕地 识别 的 “窗口 期 " ;从 整体 来 
看 ,中 低 覆 盖 度 葡萄 林 与 非 植 被 ,由 于 区 域内 植被 覆 
盖 度 低 ,NDVI 值 整体 较 低 ,葡萄 林 相 对 非 植被 和 稀 
PA RE ee eg , ,加 之 受 人 类 活动 ( 打 梢 ) 的 影响 
较 大 ,使 其 季节 变化 明显 , 即 中 低 覆 盖 度 葡萄 林 的 
NDVI 差 值 速率 变化 较 裸 地 和 稀 玖 草地 明显 。 因 
此 ,充分 结合 葡萄 林 时 序曲 线 中 NDVI 值 及 其 差 值 
速率 ,可 以 更 有 效 地 提取 中 低 履 盖 度 葡萄 林 信 息 。 

对 于 高 覆盖 度 葡萄 林 和 其 他 高 覆盖 度 植被 ， 
NDVI 都 呈现 出 生长 季 明 显 的 快速 上 升 趋势 。 小 麦 


NDVI 值 较 低 ,从 3 月 下 旬 开 始 呈 上 升 趋势 ,到 6 月 
未 达到 最 高 值 ,其 后 NDVI 值 整体 较 高 且 无 较 大 波 
动 ,在 11 月 中 旬 进 入 冬季 后 叶片 脱落 、 误 雁 ,NDVI 
值 降低 ; 玉米 在 5 月 下 旬 开 始 迅速 生长 ,7 一 9 月 
NDVI 曲线 平缓 并 保持 较 高 水 平 ,9 月 收割 后 NDVI 
值 迅 速 降低 。 从 图 3 发 现 ,高 覆盖 度 葡 萄 林 与 牧草 
的 NDVI 时 序曲 线形 状 相似 , 且 季 节 变 化 相近 ,单纯 
利用 NDVI 值 及 其 差 值 变化 速率 很 难 将 其 分 开 。 总 
体 看 来 ,高 覆盖 度 葡 萄 林 与 其 他 高 覆盖 度 植被 的 
NDVI 曲线 有 多 处 相交 , 差 值 速 率 较为 相近 ,然而 
NDVI 时间 序列 曲线 的 积分 值 能 有 效 放 大 地 物 类 型 
的 特征 ,显著 提升 葡萄 林 与 其 他 地 物 的 差异 并 可 以 
减少 影像 中 “ 同 谱 异 物 ” 和 “ 同 物 异 谱 ” 的 影响 。 
此 ,可 利用 相 邻 时 相 的 NDVI 曲线 积分 值 将 高 覆盖 
度 葡萄 林 进 行 识别 。 
3.2 J-M 距离 及 最 佳 时序 NDVI AA 

依据 研究 区 葡萄 林 生 长 的 物候 差异 及 其 完整 生 
长 周期 的 NDVI 曲线 特征 ,可 获得 不 同时 序 NDVI 的 
组 合 数据 ,通过 计算 不 同时 序 NDVI 组 合 数据 的 J-M 
距离 ,确定 葡萄 林 提 取 的 最 佳 时 序 ,减少 数据 的 元 余 


的 季节 变化 十 分 明显 ,呈现 出 双 峰 曲线 特征 ,其 原因 
是 7 月 小 麦收 割 后 会 有 蔬 荣 种 植 ; 林 地 在 休眠 期 内 


量 , 以 提高 提取 精度 。 
JM 上 距离 计算 结果 ( 表 5) 表 明 ,和 葡萄 林 与 其 他 


lp -一 中 低 覆 盖 度 葡萄 林 一 高 覆盖 度 葡萄 林 一 一 草地 -o- 小 麦 一 一 牧草 一 -林地 -x- 玉米 
0.8 F Cm ——— x ———X 
0.6 
$ 
z 0.4 F 
2 
0.2 F 
一 一 一 一 一 X ay mn 
0 e. 1 1 1 1 1 j 
03-26 04-28 05-26 06-16 07-06 07-30 08-16 09-02 10-12 11-14 12-05 
-0.2 日 期 /月 -日 
图 3 NDVI 时 间 序 列 变化 曲线 
Fig.3 Variation of time series of NDVI 
X5 Jeffreys-Matusita 距离 计算 结果 
Tab.5 Calculated results of Jeffreys-Matusita distance 
NDVI 数据 编号 组 合 GPT-YM GPT- MC&SC  GPT- XM GPT-FHL  ZDPT-GPT  ZDPT-SCD ZDPT - LD 
2-3-4 1.980 1.990 1.978 1.996 1.69 1.855 1.838 
2-3-4-5 2.000 1.999 1.998 1.998 1. 790 1.965 1.868 
2-3-4-5-6-7-8 1.999 2.000 1.998 2.000 1.988 1.962 1.984 
2-3-4-5-6-7-8-9 2.000 2.000 2.000 2.000 2.000 2.000 2.000 
1-2-3-4-5-6-7-8-9 2.000 2.000 2.000 2.000 2.000 2.000 2.000 
2-3-4-5-6-7-8-9-10 2.000 2.000 2.000 2.000 2.000 2.000 2.000 
2-3-4-5-6-7-8-9-10-11 2.000 2.000 2.000 2.000 2.000 2.000 2.000 


IRRE ;LD 表示 裸 地 。 


注 :GPT 表示 高 覆盖 度 和 葡萄 林 ;ZDPT 表示 中 低 覆 盖 度 葡萄 林 ;YM 表示 玉米 ;MC&SC 表示 牧草 和 蔬菜 ;XM 表示 小 麦 ;FHL 表示 防护 林 ;SCD K 
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地 物 的 分 离 性 随时 序数 据 的 增加 而 增 大 ,但 当时 序 
数据 组 合 为 4 一 10 月 时 ,J-M 距离 值 为 2, 可 分 离 性 
已 达到 最 住 水 平 。 因 此 ,根据 J-M 距离 计算 结果 , 确 
定 NDVI 数据 组 合 4 一 10 月 的 8 个 时 相 为 最 佳 时 序 
数据 组 合 ,并 将 此 时 序 用 于 后 面 的 波段 运算 ,计算 天 
值 与 S$ 值 。 
3.3 葡萄 林 识 别 模型 构建 

不 同 地 类 的 归 一 化 指数 分 析 发 现 :(D) NDVI 指 
数 可 以 将 植被 覆盖 地 类 分 为 高 覆盖 植被 类 和 低 履 盖 
植被 类 及 裸 地 ,水体 ;@) 低 覆 盖 植 被 类 中 ,中 低 覆 盖 
度 葡萄 林 的 NDVI 指数 均值 最 高 ,而 且 季节 变化 明 
显 ,可 采用 NDVI 指数 结合 差 值 速率 (K) 区 分 中 低 
窗 盖 度 葡萄 林 与 其 他 2 种 地 类 ;@) 高 覆盖 植被 类 
中 ,高 覆盖 度 葡萄 林 与 耕地 、 林 地 、 牧 草 覆 盖 度 相对 
较 高 ,其 生长 旺季 的 NDVI 值 过 大 易 饱 和 ,利用 相 邻 
时 相 的 面积 曲线 值 放 大 各 地 类 的 特征 值 ,可 达到 区 
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NDVI ,<0 


0730 d 
Kos oses » -0.001 
NDVI pı > 0.075 


NDVI,,,, < 0.1 


0428 ^ 


» 0.461 


分 地 类 的 目的 。 

结合 NDVI 曲线 分 析 和 数据 挖掘 ,首先 ,可 用 植 
被 指数 NDV Loo: 进行 闵 值 分 割 为 NDVI > 0. 461 
(耕地 、 林 地 牧草 、 高 覆盖 度 和 葡萄 林 ) 和 NDVlow, < 
0. 461 ( 中 低 窗 盖 度 葡萄 林 、 草 地 、 裸 地 水 体 ) ;其 
次 ,运用 Clementine 中 C5. 0 决策 树 分 类 算法 分 别 建 
立 中 低 覆 盖 度 葡萄 林 和 高 覆盖 度 和 葡萄 林 的 决策 树 ， 
并 获取 规则 ;最 后 ,结合 葡萄 林 生 长 物候 期 和 反复 试 
验 ,判断 NDVIows 为 葡萄 提取 的 关键 时 期 之 一 。 由 
于 该 时 期 葡萄 林 刚 从 地 里 挖 出 ,与 裸 地 和 政 草 地 的 
光谱 信息 十 分 相似 。 因 此 ,在 葡萄 提取 过 程 中 可 加 
入 这 一 关键 时 期 ,减少 葡萄 小 苗 与 裸 地 和 和 草地 的 混 
消 。 最 终 , 经 过 反复 试验 ,得 到 葡萄 林 提 取 规 则 
(图 4): 

中 低 覆 盖 度 葡萄 林 :NDVIwu <0. 461 ,NDVIS; >0， 
Kops oos > —0. 001, NDVI,, 20.075, NDVI,4, «0.1; 
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0526-0616 


» 9.518 


其 他 (小 麦 . 林地、 牧草 ) 


S 


上 0706-0730 


< 16.307 


| 其他 


| rcm R NDVI, < 0.07 m 


| 其 他 


高 覆盖 度 葡萄 林 


图 4 决策 树 分 类 流程 及 参数 


Fig.4 Classification process and parameters of the decision tree 


T3181 si HER] ARP : NDV Ig, > 0. 461, Sosz oss S 
9.518, Soo .ono € 16. 307, Sq. > 4. 746, 
Sono osi > 11. 309, NDVIg4, «0.07, 

3.4 分 类 结果 与 精度 验证 

采用 上 述 方法 对 红寺堡 区 2017 年 的 葡萄 种 植 
区 进行 提取 ,获得 的 分 类 结果 如 图 5 所 示 。 考 虑 到 
本 文 的 研究 对 象 为 葡萄 林 , 因 此 将 其 他 地 物 类 型 
合并 。 

由 图 5 可 以 看 出 ,研究 区 内 葡萄 种 植 区 多 为 不 
连续 的 小 面积 区 块 , 较 大 面积 区 块 较 小 ,破碎 区 块 较 
多 。 且 葡萄 种 植 区 集中 在 研究 区 中 部 , 罗 山 脚下 片 
区 分 布 面积 较 大 且 较 为 连贯 。 其 次 在 东部 和 南部 有 
较 少 分 布 。 在 缺少 2017 年 葡萄 种 植 面积 统计 数据 
的 情况 下 , 本 文 将 引用 最 新 数据 提 及 的 种 植 面 积 进 
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图 5 葡萄 林 提 取 结 果 


Fig.5 Extracted results of grape 


表 6 决策 树 分 类 精度 评价 
Tab.6 Evaluation on classification accuracy of the 


decision tree 


类 别 PU/% UA/% 
中 低 覆 盖 度 葡萄 林 92.10 90. 32 
高 覆盖 度 葡萄 林 91.20 90.21 
其 他 植被 95.06 97.58 
非 植 被 92.92 87.50 
总 精度 94. 49 
Kappa 系数 0.90 


注 :PU UA 分 别 代表 各 类 的 生产 精度 和 用 户 精度 。 


行 比较 。 有 关 数 据 统计 ,截至 2013 年 底 , 红 寺 堡 区 
葡萄 累计 种 植 总 面积 达到 84 km? P 。 分 类 结果 显 
示 , 分 类 提取 葡萄 林 总 面积 为 71.3 km, FREK 
度 葡 萄 林 面 积 为 56 km ;高 覆盖 度 葡萄 林 面 积 ; 
14.7 km 。 提 取 相 对 误差 为 15.1% 。 

将 人 机 交互 决策 树 生 成 的 分 类 结果 同 野外 采样 
验证 样本 集 数据 作对 照 ,建立 混淆 矩阵 ,选取 生产 者 
精度 ,用户 精度 Kappa 系数 和 总 精度 4 个 指标 对 分 
类 结果 进行 验证 与 评价 。 为 了 对 比 决 策 树 模型 构建 
的 优 劣 程度 ,本 研究 选择 同样 训练 样本 ,采用 支持 向 
量 机 分 类 法 进行 对 比 , 结 果 如 表 6 所 示 。 结 果 表 明 ， 
基于 决策 树 分 类 的 总 体 精度 为 94.49% ,Kappa 系数 
为 0. 90。 其 中 ,中 低 覆 盖 度 葡萄 林 生 产 精 度 为 
92. 10% ,用 户 精 度 为 90. 32% ; 高 覆盖 度 葡 萄 林 生 
产 精 度 为 91.20% ,用 户 精 度 为 90.21% 。 


4 讨论 与 结论 


葡萄 生长 季 与 其 他 植被 较为 近似 ,使 用 某 一 时 
相 的 遥感 影像 提取 葡萄 种 植 区 域 , 极 易 发 生 混淆 ,大 
幅度 降低 提取 精度 。 而 选择 最 佳 NDVI 时 间 序 列 曲 
线 特征 并 结合 葡萄 林 生 长 物候 期 和 种 植 规律 的 分 类 
识别 方法 ,对 葡萄 林 的 种 植 信 息 提 取 精 度 较 高 ,结果 
满足 需求 。 覆 盖 葡 萄 林 生 育 期 的 时 间 序 列 遥 感 数据 
完整 地 反映 了 和 葡萄 林 长 势 随时 间 的 变化 ,利用 恰当 
的 数学 方法 表现 特殊 的 物候 特征 ,能 够 有 效 将 其 与 
其 他 地 物 区 分 开 ,识别 出 葡萄 林 种 植 地 块 像 元 。 该 
方法 从 地 物 的 NDVI 曲线 的 变化 特征 入 手 ,根据 不 
同 履 善 度 的 葡萄 林 换 索 NDVI 曲线 的 数学 模型 放大 
物候 特征 ,以 数据 挖掘 的 方法 寻找 特殊 节点 ,该 方法 
所 需 参 数 较 少 ,自动 化 程度 高 ,适用 性 强 ,适宜 向 更 
大 范围 的 葡萄 林 种 植 区 域 推广 ,有 和 较 高 的 实用 价值 。 
总 而 言 之 ,具有 较 高 时 间 和 空间 分 辩 率 优势 的 GF - 


1/WFV 遥感 数据 能 够 较 好 地 重建 葡萄 林 生 长 季 
NDVI 曲线 ,反映 其 变化 特征 ,因此 ,通过 选择 葡萄 
林 生 长 期 最 佳 时 相 的 遥感 数据 ,然后 利用 时 序 植被 
指数 进行 葡萄 林 分 类 ,可 以 有 效 减 小 “异物 同 谱 ? 现 
象 的 干扰 。 但 从 分 类 结果 看 ,提取 葡萄 林 图 斑 完 整 
性 较 低 ,其 主要 原因 为 红寺堡 区 不 同 的 地 貌 .土壤 条 
件 等 ,导致 葡萄 林 生 长 速度 .长 势 状态 等 都 有 一 定 的 
差异 ,影响 了 葡萄 林 提 取 的 准确 性 。 因 此 ,在 未 来 的 
研究 中 ,应 采用 高 分 辩 率 影像 融合 的 方法 ,结合 使 用 
面向 对 象 分 类 法 ,有 针对 性 的 设置 参数 ,将 葡萄 林地 
块 进行 划分 。 

本 文 研 究 得 出 :充分 利用 葡萄 林 的 生长 物候 期 
和 种 植 规律 的 NDVI 时 序 特 征 , 针 对 不 同 覆盖 度 的 
葡萄 林 选 择 恰 当 的 分 类 方法 ,可 准确 有 效 提取 葡萄 
林 。 基 于 2017 Æ 4—9 H GF - 1/WEV 时 间 序 列 数 
据 并 结合 实地 采集 数据 信息 ,利用 阔 值 分 割 方法 将 
研究 区 内 的 高 、 低 覆盖 度 植 被 进行 划分 ,通过 进一步 
控 气 与 分 析 葡 萄 的 提取 窗口 期 ,发现 中 低 覆 盖 度 葡 
萄 林 可 利用 NDVI 及 其 差 值 速率 进行 有 效 识别 ;高 
窗 羡 度 葡萄 林 可 以 利用 NDVI 及 其 相 邻 时 相 曲 线 面 
积分 别 进行 有 效 识别 。 

运用 J-M 距离 判别 方法 ,研究 了 和 葡萄 林 与 其 他 
地 物 类 型 区 分 性 能 最 佳 的 波段 组 合 , 即 得 到 了 葡萄 
林 提 取 的 最 佳 时 间 序 列 。 研 究 指 出 :4 一 9 月 的 时 间 
序列 ,由 4 月 28 日 5 月 26 日 6 月 16 日 7 月 2 日 、 
7 月 30 日 8 月 16 日 ,9 月 2 日 共 7 个 时 相 组 成 葡萄 
林 提 取 的 最 佳 时 相 , 使 得 葡萄 林 和 其 他 地 类 的 区 分 
性 较 好 且 数据 宛 余 降 低 。 

基于 特征 波段 的 决策 树 分 类 方法 提取 和 葡萄 林 ， 
深入 分 析 和 挖掘 光谱 信息 ,有 效 提 取 了 研究 区 葡萄 
分 布 且 达到 针对 性 强 、 分 类 精度 明显 提高 的 目的 。 
对 识别 模型 及 提取 精度 进行 了 验证 评价 ,提取 总 精 
度 达 到 了 94% ,从 而 证 明了 本 文 葡萄 林 识 别 模 型 的 
有 效 性 和 可 靠 性 。 

研究 中 存在 的 不 足 之 处 ,主要 表现 为 :对 于 空间 
分 辩 率 问题 导致 的 混合 像 元 没有 做 处 理 , 在 一 定 程 
度 上 影响 了 提取 误差 ,在 葡萄 林 提 取 中 ,部 分 的 道路 
像 元 被 划分 到 葡萄 林 像 元 中 ;分 类 结果 显示 ,研究 区 
内 部 分 乡间 道路 旁 的 低 矮 灌 从 ,与 葡萄 种 植 地 区 的 
时 间 序 列 光谱 特征 十 分 类 似 , 造 成 混淆 。 考 虑 到 道 
路 在 整个 研究 区 内 所 占 面积 较 小 , 且 该 误差 对 葡萄 
种 植 面积 提取 的 影响 不 大 ,所 以 ,在 误差 允许 的 范围 
内 ,未 来 研究 中 将 予以 改进 。 
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Grape Recognition Model Based on GF — 1/WFV Time Series: 
A Case Study in Hongsibu District of Ningxia 
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Abstract: Grape is one of the most widely distributed fruit tree species, and its accurate spatial distribution is of 
great significance for the management and development of grape planting and wine industry. In this study ,the Hong- 
sibao in Ningxia was taken as the study area to obtain the time series of spectrum data in grape growing season based 
on the Gaofen-1 satellite (GF — 1L/WFV) ,the J-M distance analysis was used to analyze the normalized difference 
vegetation index ( NDVI) of the grape plots, and the best recognition phase was determined. The NDVI values of the 
best phase , adjacent to the value difference rate and curvilinear integral training samples were put into the Clementine 
data mining software. The extract grape model of decision tree was developed by using C5. 0 decision tree classification 
algorithm and combining with the expert experience method. The results suggested that the established recognition 
model could meet the needs of grape recognition , but the accuracy was different from different grape plots due to their 
different coverage. The overall accuracy based on the decision tree classification was 93. 7196 ,and the Kappa coeffi- 
cient was 0. 91. In which the production precision and users' precision of the grape plots with moderate or low coverage 
were 90. 82% and 88.5646 ,and those with high coverage were 92. 44% and 91. 1846 ,respectively. 
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